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INTRODUCTION 

Since the discovery of the High-T c superconductors, [1], the Hubbard model has 
been the subject of intense investigations following Anderson's proposal [2] that the 
model should capture the essential physics of the cuprate superconductors. From the 
earlier attempts to obtain the magnetic phase diagram on the square lattice (for an 
overview see the book by Mattis [3]) one can deduce that antiferromagnetic order 
exists in the vicinity of the half-filled band whereas ferromagnetic ordering might take 
place in the phase diagram for strong repulsive interaction strength and moderate hole 
doping of the half-filled band. Obviously antiferromagnetic and ferromagnetic orders 
compete in this part of the phase diagram. More recent calculations [4] established 
that the ground state of the Hubbard model on the square lattice shows long-ranged 
antiferromagnetic ordering with a charge transfer gap. However, the problem of mobile 
holes in an antiferromagnetic background remains mostly unsolved. Suggestions for a 
very wide ferromagnetic domain in the phase diagram based on the restricted Hartree- 
Fock Approximation have been made by several authors [5] on the cubic lattice, and 
on the square lattice [6-8]. This domain appears for large interaction and moderate 
hole doping in which case the Hartree-Fock Approximation ceases to be controlled. 
Within this framework one expects to obtain reliable results for moderate U where the 
paramagnetic phase is indeed unstable towards an incommensurate spin structure at a 
critical density n c (U) [9]. The Gutzwiller Approximation (GA) [10-12] has been applied 
[13], even for large U, yielding results similar to the Hartree-Fock Approximation. 
However, for large U, a ferromagnetic domain appears only if the density is larger 
than some critical value. In the Kotliar and Ruckenstein slave boson technique [14] 
the GA appears as a saddle-point approximation of this held theoretical representation 
of the Hubbard model. In the latter a metal-insulator transition occurs at half-filling 
as recently discussed by Lavagna [15]. The contribution of the thermal fluctuations 
has been calculated [16] and turned out to be incomplete as this representation, even 
though exact, is not manifestly spin-rotation invariant. Spin-rotation invariant [17] and 
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spin and charge-rotation invariant [18] formulations have been proposed and the first 
one was used to calculate correlation functions [19] and spin fluctuation contributions 
to the specific heat [20]. Comparisons of ground state energy with Quantum Monte- 
Carlo simulations, including antiferromagnetic ordering [21] and spiral states [22], or 
with exact diagonalisation data [23] have been done and yield excellent agreement, and 
a magnetic phase diagram has been proposed [24]. 

Such calculations are not restricted to the square lattice and can be performed, for 
example, on the honeycomb lattice, where a (semi-) metal to insulator transition takes 
place at half-filling when the interaction strength is increased [25]. In this paper we first 
compare the slave-boson estimate of the critical interaction strength with the Quantum 
Monte-Carlo simulations of Sorella and Tossati [25]. In a second part we determine the 
phase boundary of the paramagnetic state on the square lattice by determining the 
instability line of the latter with respect to an arbitrary long-ranged incommensurate 
spin order at finite temperature, thus extending the calculations we performed at zero 
temperature [26]. 



HONEYCOMB LATTICE RESULTS 



Considering the honeycomb (HC) lattice instead of the square lattice allows for 
a further study of the interaction induced Metal-to- Insulator Transition (MIT). It has 
been shown by Quantum Monte-Carlo simulations [25] that, at half-filling, the system 
undergoes a MIT as the interaction strength is raised up. This follows from the fact 
that we meet here a semi-metal in which, as far as we consider the Stoner criterion, no 
long-ranged magnetic order is expected for small U . On the other hand the HC lattice 
is bipartite and can obviously support un-frustrated anti-ferromagnetism at large U . 
The two considerations thus yield a scenario for the MIT at half-filling. By raising up 
the interaction long-ranged anti-ferromagnetic order will set in and a gap will open at 
a critical interaction strength, thus yielding a transition from a metallic to an insulator 
behavior. Another possibility is provided by the Brinkman-Rice mechanism. Here 
raising up the interaction leads to a diverging effective mass [11, 12] at a critical value 
of U, and a gap opens in the charge excitation spectrum [15, 27, 28] at the Fermi energy. 
It thus describes a Mott insulator. These two scenarii can be addressed in the slave 
boson mean-field approach. 

On the HC lattice, the free electron dispersion relation is given by: 

t% = 2tv | 1 + exp (— ik ■ bi) + exp (— ik ■ b 2 ) | v = ±1 (1) 

where, a being the lattice constant, b\ = (3, y/3)a/2, b 2 = (3, — y/3)a/2 and the pseudo- 
fermions acquires the effective dispersion relation E^ au = z^t^ + f3 — fi where Zq 
denotes the inverse effective mass of the quasi-particles: 

*l = ^f( e + J ) 2 ! ? ■ ( 2 ) 

The paramagnetic mean-field energy reads: 

F = -T J2 In (1 + e- (E ^ IT) ) + Ud 2 + a(e 2 + p 2 + d 2 - 1) - /3 (2d 2 + p 2 ) (3) 

where e 2 ,j>Q,G? 2 denote respectively the average occupation number of empty, singly- 
occupied, and doubly-occupied sites ; a and f3 are the Lagrange multipliers enforcing 
the constraints [18]. When staggered order is allowed, the effective dispersion relation 
for the pseudo-fermions becomes: 
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*W = A) " + a^^T^^T a = ±1 (4) 



where 

z _ 1 Po(e + d) + (Jp(e - J) 1 

\Jl - d 2 - V2 \Jl~e 2 - iB^sE 

and the free energy reads: 

F = -T J2 W 1 + e" ( ^ /T) ) + ^ 2 + a(e 2 +^ + p 2 + J 2 -l)- / 3 (2J 2 +^ +p 2 )_2p /3-p 

(6) 

We minimize both free energies eq. (3), eq. (6) at zero temperature. The resulting 
free energies are displayed in fig. fa. The agreement with the Quantum Monte-Carlo 
values of Sorella and Tossati [25] is excellent as long as no magnetic order occurs, even 
though no quantum fluctuations are included in the calculation. The differences range 
from 5% at U = it to 11% at U = 8t. 




Figure 1. a) Ground state energy of the Hubbard Model on the honeycomb lattice at half- 
filling as a function of the interaction strength in the paramagnetic state (Curve A) and in 
the anti-ferromagnetic state (curve B). 

b) Density of doubly occupied sites for the Hubbard Model on the honeycomb lattice at 
half-filling as a function of the interaction strength in the paramagnetic state (Curve A) and 
in the anti-ferromagnetic state (curve B). 



In fig. lb we show the double occupancy in both paramagnetic and anti- 
ferromagnetic phases. In the latter the double occupancy is seen to be reduced from 
its value in the paramagnetic phase in the intermediate coupling regime. In the strong 
coupling regime, i. e. above the Brinkman-Rice point, it will be the opposite and 
d 2 will be proportional t/U. As compared to the simulations [25], our results for the 
double-occupancy are found to be slightly smaller, especially for the largest available 
values of U . In fig. 2a we locate the metal to insulator transition point by showing the 
staggered magnetization m as a function of U . It turns out that m starts to deviate 
from at U ~ 3. It, a value that is fairly smaller than the one obtained in Quantum 
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Monte-Carlo simulations (4.5 ±0.5) but still reasonably close to it. It represents a size- 
able improvement of the Hartree-Fock value (2.2). It also rules out the Brinkman-Rice 
scenario, as the critical U is 12. 6t. The magnetic gap is given by twice the Lagrange 
multiplier f3. It is displayed in fig. 2b. It starts to deviate from at about U/t ~ 3.f 
and grows fairly linearly with U — U c in the vicinity of the transition and is proportional 
to U in the strong coupling limit. 




Figure 2. a) Square of the staggered magnetization of the Hubbard Model on the honeycomb 
lattice at half-filling as a function of the interaction strength. 

b) Magnetic gap of the Hubbard Model on the honeycomb lattice at half-filling as a function 
of the interaction strength. 



MAGNETIC PHASE DIAGRAM 



In the spin-rotation invariant slave boson mean-held theory of the Hubbard model 
on the square lattice the ground state mean held energy is given by: 



F = -T V In (1 + e- {E ^' T) ) + Ud 2 + a(e 2 + p 2 + d 2 - 1) - f3 (2d 2 + p 2 ) (7) 



k,<T 

with the pseudo-fermion dispersion relation E% = z^t^ + f3 — fi . Minimizing F with 
respect to e, p 0} d } a and f3 delivers the saddle-point. If, using the constraints, one 
expresses e and p as functions of d and the density, one immediately obtains the free 
energy (7) as expressed in the Gutzwiller approximation. 

As discussed by many authors, those saddle-point conditions can be cast into a 
single one: 

U l-(e + d) 2 

u c = (e + dy-p {e + d) (8) 

where 

/oo -/2 
deeN(e)f F (^(e- Mf )) (9) 
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where N(e) denotes the density of states per spin, and /i e jj = £s_M_ As usua l e an( J J 

are related by e 2 — d 2 = 8 , 8 being the hole doping. We here adopt a new attitude in 
front of the saddle-point condition eq. (8). Instead of solving directly eq. (8), we hrst 
choose a value for e and the density, and we determine numerically fi e ff- This being 
settled, we calculate U c via its definition, and we finally determine U with the help of 
eq. (8). The result of this investigation is displayed in fig. 3. We show z as a function 
of U at the finite temperature | for several densities. 




Figure 3. Solution of eq. (8) for several densities. They run on the right border from 8 = 5% 
for the upper curve down to 8 = for the lower curve by steps of 0.5%. 

Contrary to the zero temperature case, where our approach does not yield results 
differing from the standard one, we see that, for very small hole doping, we obtain 
several solutions. Moreover we obtain a line for a hrst order phase transition. The 
latter is fairly linear and is ranging between < 8 < 0.03 and 8.5 < U/t < 10. In 
detail the latter is obtained by selecting the value of z which is leading to the lowest free 
energy which, at half-filling, takes the value —Tin 4. This also gives a rough estimate 
of the Neel temperature. At half-filling and for large U, the mean-held free energy in 
the anti-ferromagnetic state is: 

F AF = -\j y (f0) 

Thus long-ranged anti-ferromagnetic orders disappears for 

T > (it 2 ) /(U In 4). (ff) 

However this approach only brings a hrst order transition, the result (fO) being essen- 
tially insensitive to the temperature. Moreover it misses the spin degrees of freedom, 
as In (2) should appear in eq. (f f ) instead of In (4). However the above-mentioned dif- 
ficulties seem to be restricted to the close vicinity of half-filling. Indeed this approach 
is expected to make sense when the temperature is much smaller than the effective 
band width, which is not the case when z is vanishingly small. It thus deserves a 
further study of the stability of the paramagnetic state with respect to an incommen- 
surate long-ranged magnetically ordered state in the overall phase diagram. Moreover 
Moreo recently [29] calculated the uniform spin susceptibility xs in a Quantum Monte- 
Carlo simulation. She obtained that, for moderate interaction strength (U/t = 4) xs 
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uniformly decreases as a function of the hole doping, whereas for strong interaction 
(U/t = 10) it hrst passes through a maximum, which is also seen in exact diagonal- 
ization for the t — J model at J/t = 0.4 [30], a result that can not be obtained in 
RPA calculations. There would then be a crossover between small and strong coupling 
behavior, even though the crossover line is not yet established [31]. Let us than see 
what our approach can say about it. 

Considering the quantum fluctuations around the saddle-point eq. (7) brings out 
2 distinct channels, spin symmetric and spin antisymmetric, providing a microscopical 
basis for a Landau- Fermi liquid theory. It is by considering the antisymmetric channel 
that one obtains the spin susceptibility. In the static limit the divergence of the latter 
yields the phase boundary of the paramagnetic state and allows to give a magnetic 
phase diagram of the Hubbard model on the square lattice at finite temperature. As 
derived by Li, Sun and Wolfle [19], the spin susceptibility is given by: 



Xs{q,u 



Xo(<?» 



l + Arxo(<?» + Aixi(?,w) + A 2 (xi(q,u) - xo(q,u)x2(q,u)) 



12) 



where: 



Xn(q,u) = - J2( t p + t p+Q) nG o{p,i^n)G {p + q,oj + iu n ) n = 0,l,2 (13) 



p,W„,<T 



and the others symbols are defined in ref. [19]. Setting u = we obtain the instability 
line of the paramagnetic phase towards an incommensurate magnetically ordered state. 
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Figure 4. a) Static Spin Susceptibility as a function of the hole doping for U/t = 4 (curve 
A) and U/t = 25 (curve B) at temperature t/6. The thick lines indicate that the calculation 
has been performed in a domain where the homogeneous paramagnetic state is stable, 
b) Phase diagram of the Hubbard Model on the square lattice at temperature t/6 as obtained 
from our approximation. 



At very small temperature (t/T > 10) the instability line is essentially unchanged 
compared to the zero temperature result [26]. But there is a sudden change at about 
t/T = 8 and we display the result at t/T = 6. The very wide magnetic domain reduces 
to a small pocket around half-filling and U/t = 10, and otherwise the paramagnetic 
state remains stable. This finds an explanation in term of the magnetic gap that 
appears to be very small in the spiral state at zero temperature as discussed in [23, 24] 
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in the largest part of the phase diagram. Allowing for a small temperature mixes the 
magnetic sub-bands and as a result the long-ranged magnetic order disappears. This 
is again a result that can not be obtained in RPA calculations, as there one obtains a 
very large magnetic gap that does not allow for a mixing of the 2 magnetic sub-bands 
at low temperature. 

In the thus stable paramagnetic state one can then compute the uniform magnetic 
susceptibility Xs- The results are displayed in fig. 4a. In the limit of small U the 
above expression for the magnetic susceptibility eq. (12) reduces to the RPA result 
and we obtain that Xs is monotonically decaying upon doping. However for large U 
Xs first passes over a maximum before decaying as well. This maximum is reminiscent 
from the one that is leading to ferromagnetic instability at zero temperature [24]. The 
obtained behavior is in agreement with the numerical simulation of Moreo [29]. Our 
result implies that it is not necessary to invoke next-nearest neighbor hopping in order 
to obtain such a behavior. That our result for Xs yields too large numbers is not 
very surprising as we are calculating in the vicinity of a magnetic instability and a 
better calculation is required in order to obtain a better agreement. All together our 
results lead to the phase diagram that is displayed in fig. 4b. It consists of 3 regions. 
The largest part, including half-filling for U/t < 3, and all densities for U/t > 21, is 
paramagnetic. This means that allowing for small temperature effects is sufficient to 
suppress the magnetic instability that was present at T = [26]. It only remains at 
intermediate coupling and small doping. In contrast to the T = case, a new phase 
separated domain appears for strong coupling and small doping. It is characterized by 
a negative compressibility and is paramagnetic. The 2 components entering a Maxwell 
construction would be paramagnetic with different densities and effective masses, one 
of them being infinite. 

CONCLUSION 

In this paper we used the Kotliar-Ruckenstein slave boson technique to treat the 
strong correlation of the Hubbard model. In a first part we discussed the metal to 
insulator transition that is occurring on the honeycomb lattice at half-filling. We find 
the critical value for the interaction strength to be U/t = 3.1. Above it a gap opens 
in the spectrum which behaves like ill — U c ) in the vicinity of the transition, and like 
U for strong coupling. In a second part we calculated the phase diagram of the Hub- 
bard model on the square lattice at finite temperature. The Mermin- Wagner theorem 
is fulfilled in the largest part of the phase diagram. In the strong coupling regime, 
the uniform magnetic susceptibility shows a maximum in its doping dependence in 
agreement with numerical simulations. 
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